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ABSTRACT 

This is the third in a series of papers reporting on a large reverberation-mapping campaign aimed to study the 
properties of active galactic nuclei (AGNs) with high accretion rates. We present new results on the variability 
of the optical Fe II emission lines in 10 AGNs observed by the Yunnan Observatory 2.4m telescope during 
2012-2013. We detect statistically significant time lags, relative to the AGN continuum, in nine of the sources. 
This accurate measurement is achieved by using a sophisticated spectral fitting scheme that allows for apparent 
flux variations of the host galaxy, and several narrow lines, due to the changing observing conditions. Six 
of the newly detected lags are indistinguishable from the H j3 lags measured in the same sources. Two are 
significantly longer and one is slightly shorter. Combining with Fe II lags reported in previous studies, we 
find a Fe II radius-luminosity relationship similar to the one for H/3, although our sample by itself shows no 
clear correlation. The results support the idea that Fe II emission lines originate in photoionized gas which, for 
the majority of the newly reported objects, is indistinguishable from the H/3-emitting gas. We also present a 
tentative correlation between the lag and intensity of Fe II and H/3 and comment on its possible origin. 

Subject headings: galaxies: active — galaxies: nuclei — galaxies: Seyfert — methods: data analysis — 
quasars: emission lines 


1. INTRODUCTION 

Most active galactic nuclei (AGNs) show prominent Fe II 
emission lines in their spectra. The lines appear in several 
broad bands that represent thousands of individual transitions. 
The stronger bands cover the ranges of 4000-5400 A (here¬ 
after optical Fe II lines), 28 00-3500 A and 2000-2600 A 
(hereafter UV Fe II lin es; e.g. JWills et alJl985tlSulentic et all 
120001: iHu et all 120081 and references therein). The strongest 
Fe II lines (relative to H/3) are observed in narrow-line 
Seyfert 1 galaxies (NLSls, see, e.g., lOsterbrock & Poggej 
1 1 985b[Boiler et alJll996l) . Such objects are characterized by: 
1) narrow (< 2000 km s -1 ) broad emission lines; 2) weak [O 
III] line; 3) and steep hard X-ray spectrum. These objects are 
found usually a t the extreme end of the so-called Eigenvec¬ 
tor 1 sequence jBoroson & Greenlll992tlSulentic et aUll2000fc 
IShen & Holl2014l) indicating high Eddington ratios and sev¬ 
eral other properties that are not fully understood. 

Some NLSls show evidence for super-Eddington accre¬ 
tion, with Lboi/TEdd > E We refer to these objects as super- 
Eddington ac creting massive black h oles (SEAMBHs). Our 
earlier work (IWang et al.ll2013L 120141) show that such objects 
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are potentially a new kind of standard candle for cosmology. 
To test this idea, and to study in more detail the physical 
properties of extreme NLSls, we have initiated a large rever¬ 
beration mapping (RM) campaign to measure, accurately, the 
black hole (BH) mass of such sources. The initial results of 
this campaign hav e been publ i shed i n iDu et alJ (12014 here¬ 
after Paper I) and I Wang et all (12014 Paper II). The present 
paper, the third in this series, is dedicated to the study of Fe 
II emission lines in our sample of SEAMBH candidates. A 
fourth paper, submitted in parallel to this one, presents new 
data on H/3 time lags and BH mass in five extreme SEAMBHs 
and focuses on the modification of the Rm^-L relation ship in 
AGNs (e.g. JKaspi et~d]l2000ll200aiBentz et allMl in the 
presence of such sources. 

The excitation mechanism of the Fe II emission in 
AGNs has been discussed in numerous publications (e.g. , 
Collin-Souffrin et al J[l980t iNetzer & Willsll 1 9831: l.lolv II 1987k 


Sigut & Pradhanl Il998t iBaldwin et al. 1 12004 iFerland et all 

2009). Most of these studies suggest an origin in the broad¬ 
line region (BLR) but the line intensities calculated so far are 
in poor agreement with most observations. Other suggestions 
connect these lines to the outer part of the central accretion 
disk (e.g.. lJolv~lll987l) but the agreement with observations 
is still poor. In particular, none of the existing models can 
explain, satisfactorily, the relative intensity of the UV and op¬ 
tical Fe II lines and the observed Fe II spectrum in the range 
2000-2600 A. Line variability is an important tool in such 
studies since it can indicate, given a measured time lag rel¬ 
ative to the continuum variations, the location of the Fe II- 
emitting gas. Such variations have been detected in a_num- 
ber of sources (e.g.,f Boksenbei^_&NgtzeijLl973;lMaoz_gtalJ 


19931 fK ollatsch nv et al. 2000t Vestergaard & Petersonl 2005 

Wang et all 120051 : iKuehnetalJ 120081 : IShapovalova et all 
20121) , but robust lag measureme nts were not obtai ned. This 
situation has changed, recently . iBian et al.1 (1201 Oh revisited 
the data of PG 1700+518 from Kaspi et al. (2000) and mea- 
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sured the light curve of the optical Fe II lines. They ob¬ 
tained a significant time lag relative to the 5100 A contin- 
uu m albeit with a ve r y larg e uncertainty. iRafter et~afl (120131) 
and lChelouche et all d2014T) ado pted the multivariate correla- 
tion function (MCF) scheme of C helouche & Zuckerl (1201 31) 
designed for photometric RM. IRafter et al.l (12013^ studied 
the NLS1 SDSS J113913.91+335551.1 using RM and found 
the Fe II time lag js consistent with the one of the H/3 
line. IChelouche & Zuckerl d2013l) measured optical Fe II light 
curves , and time lags, for three of the objects in Kasgietah 
( 2000 )’ one °f which is the object studies hv Bian et al. 
(2010). They also find somewhat less_significant Fe II time lag 
for three other sources. IChelouche & Zuckerl (120131) presents 
a tentative Fe II size-luminosity relation and suggest that the 
Fe II emission-region size is comparable to that of the H/3 
line. iBarth et ahl ( 2013 ) used spectroscopic measurements of 
optical Fe II lines and found lags that are 1.5 and 1.9 times 
longer than the corresponding H/3. All these studies confirm 
the photoionization origin of the Fe II lines. The success of 
the IBarth et ahl d2013l) campaign is due both to the detailed 
and frequent spectroscopic observations and a novel method 
they used to correct for the host galaxy contribution to the 
optical spectra. 

This paper presents the results obtained, so far, to measure 
the time lags of the optical Fe II lines (hereafter Fe II lines) in 
our SEAMBH campaign. We have attempted to detect such 
lags in 10 of the sources and were able to obtain statistically 
significant results in nine of them. In six of the newly mea¬ 
sured sources, the Fe II lag is entirely consistent with the H/3 
lag and in two others it is considerably longer. SectionOgives 
a brief review of the observations and data reduction. Section 
[^describes our spectral fitting method, with emphasis on host 
galaxy and narrow line subtraction which we find to be crucial 
to the analysis. More details are given in Appendix lAl Sec¬ 
tion |4] presents Fe II light curves and their analysis for both 
Fe II and H/3, and compares the results of the new H/3 lags 
to those presented in Paper I and Paper II. In Section 0 we 
plot the size-luminosity relation for Fe II and compare the lag 
and intensity of Fe II with H/3. The implications and some 
additional interpretation are also discussed. Section Ogives a 
summary of the new results. 

2 . OBSERVATIONS AND DATA REDUCTIONS 

The details of the SEAMBH campaign, including the obser¬ 
vations, data reduction and analysis, were presented in Paper 
I and Paper II. For completeness, we summarize the more im¬ 
portant points below and discuss in detail the new method of 
galaxy and narrow line subtraction. 

2.1. Sample 

Ten NLSls identified as SEAMBH candidates were ob¬ 
served, spectroscopically and photometrically, between Oc¬ 
tober 2012 and June 2013. Objects names and coordinates 
are listed in Table 1 of Paper II. H/3 time lags for three of 
the sources (Mrk 335, Mrk 142, IRAS F12397+3333) are 
presented in Paper I and for five additional sources (Mrk 
1044, Mrk 382, MCG +06-26-01 2 Mrk 486, Mrk 493) in 
Paper II. Like many other NLSls (jBoroson & Gree nl 1 19921: 
iSulentic et al .l2()00tlBorosonl2002t [Zhou et al.ll2006l e.g.). all 
the sources in our sample show strong Fe II emission lines and 
high Lboi/CEdd- This sample is, therefore, different from most 
AGNs in the local universe and none of the results presented 
below should be compared with earlier studies, like iHu et al.l 
( 120081 ) that address the general population properties. 


2.2. Spectroscopy and Data Reduction 

The spectra were obtained using the Yunnan Faint Ob¬ 
ject Spectrograph and Camera (YFOSC), mounted on the Li- 
jiang 2.4m telescope at Yunnan Observatory of the Chinese 
Academy of Sciences. A longslit with projected width of 235 
was oriented to take the spectra of the object and a nearby 
non-varying comparison star simul t aneo usly following, e.g., 
iMaoz et ah (119901) and iKasni et alJ (120001) . The comparison 
star is then used as a standard for flux calibration. For the 
Lijiang 2.4m telescope, the rotator is accurate and the track¬ 
ing is stable. Thus the object and the comparison star were 
kept within the slit during the typical 30 min. exposures. The 
distance between the object and the comparison star along the 
direction of slit width keeps less than 1 pixel (03283). 

Grism 14 was used and yielded spectra covering the wave¬ 
length range of 3800-7200 A with a dispersion of 1.8 A 
pixel” 1 . The final spectral resolution, obtained by compar¬ 
ing the width of the [O III] emission line with the one mea¬ 
sured from the Sloan Digital Sky Survey (SDSS: lYork et al.l 
120001) spectrum of the same object, is roughly 500 km s” 1 . 
All spectra are extracted in a uniform, large aperture of 835 
to minimize light losses. The flux calibration is based on the 
comparison of the object flux to the comparison star flux in 
the 2 "5 x835 apertures. However, the host galaxy flux cal¬ 
ibration is different because of the non-stellar image — the 
galaxy is extended and resolved. This can result in apparent 
flux variations due to variable seeing and mis-centering. This 
requires a different calibration procedure that has a large im¬ 
pact on the various light curves. Appendix lAl gives details of 
this procedure. 

3. LINE AND CONTINUUM MEASUREMENTS 
3.1. Fitting scheme 

The traditional method to measure the flux of the continuum 
and the broad e missio n lines, in most RM studies, is simple 
integration (e.g.. lKaspi et al.ll200Q ). A straight line is set be¬ 
tween two line-free windows to define the AGN continuum, 
and the flux of the emission line is measured by simple in¬ 
tegration above the line. This method works well for single, 
strong emission lines, e.g., H/3, because there are no (or only 
weak) other emission lines in this range and the wavelength 
window is narrow enough for the continuum to be approxi¬ 
mated by a straight line. In the case of Fe II, neither condition 
is satisfied. 

Fe II emission consists of thousands of lines that form a 
pseudo-continuum. The most prominent features in the opti¬ 
cal band are two bumps between 4500 to 5500 A: one between 
H/3 and H 7 , and the other on the red side of [O III] A5007. 
Even for quasars, whose host galaxies are extremely faint rel¬ 
ative to the AGN, it is hard to find “pure” continuum windows, 
and it is inappropriate to assume a straight line over such a 
wide wavelength range of more than 1000 A. The strong host 
galaxy contribution in most of our objects, and the contami¬ 
nation by other emission lines (e.g.. He II and coronal lines), 
make the situation much worse. Appendix |B] presents the Fe 
II light curves measured by the traditional integration method. 
All light curves have large scatter, and only less than a half 
show rough structures. A simultaneous fitting including the 
Fe II emission and all the other spectral components is neces¬ 
sary. 

Template-fitting is a widely used method to measure Fe II 
emission in single-epoch spectra of AGNs (see, e.g jHu et al.l 
120081 for a brief overview). Input-output simulations show 
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that template-fitting is a reliable measurement of Fe II emis¬ 
sion with equivalent width (EW) > 25 A in quasars dHu et alj 
120081) . It is also common to include a host galaxy compo¬ 
nent into th e fitting of Type I AGNs with stron g host contri¬ 
bution (e.g JZhou et al.ll2006h iHo & Kimll2009l) . The method 
has been adopted, recently, to measure the light curves in a 
few reverberation mapping studies, and proved to be success¬ 
ful. |Blan|etaI] (I^I^reanalyzed the spectra of PG 1700+518, 
using data from iKasni et al.l (120001). They detect ed a lag of 
20911 47 days for Fe II emission. iBarth et al l (120131) performed 
a very careful spectral decomposition including a power-law 
continuum, a host galaxy component, Fe II template, and 
other emission lines (H/3, [O III], He II, and He I). Using 
this technique they were able to obtain statistically significant 
reverberation lags of the Fe II lines in two Seyfert 1 galaxies, 
NGC 4593 and Mrk 1511. 

The spectral fitting in the present paper follows the al¬ 
gorithm of IHu et aLM2012l) . where several spectral compo¬ 
nents are fitted simultaneously by minimizing the x 2 v ' a 
Levenberg-Marquardt method. All the light curves are ob¬ 
tained directly from the results of the fitting. Before fitting, 
we correct the calibrated spectra for Galactic extinction and 
redshift. We use the Ry-dep ende nt Galactic extinc tion law 
given by ICardelli et al.l (119891) and IQ’ Donnell (11994ft and Ry 
is assumed to be 3.1. The redshift z and V-band extinction 
are taken from the_NASA/IPAC Extragalactic DatabastQ and 
ISchlaflv & Finkbeineil(l201 ll) . as listed in Papers I and II. 

The left panel of Figure Q] shows the fitting to the spectrum 
of Mrk 382 taken at Julian date (JD) 2456298, after Galac¬ 
tic extinction and redshift correction. The fitting is performed 
in the rest frame wavelength range 4150-6280 A, except two 
narrow windows around H 7 and He I A5876. The two lines 
are not blended with the major part of Fe II emission, and 
their study is beyond the scope of the present paper. We keep 
them out of the fitting to avoid introducing too many unneces¬ 
sary parameters. The observed spectra in the fitting windows 
are plotted in green, while those left out of the fitting are in 
black. The fitting includes the following components: (1) a 
single power law, (2) Fe II emission, (3) host galaxy. These 
three components, which when combined together forms the 
pseudo-continuuno are plotted in blue. (4) H/3 emission line 
plotted in magenta. (5) Broad He II A4686 emission line plot¬ 
ted in cyan. ( 6 ) Narrow emission lines plotted in orange, in¬ 
cluding [O III] AA4959, 5007, He II A4686, He I A4471, and 
several coronal lines. The summed model is plotted in red. 
The bottom panel shows the residual spectrum. Note that, al¬ 
though He I A5876 is not in the fitting window, its profile is 
well recovered after removing the Na I AA5890, 5896 (Na D) 
absorption lines from the host galaxy. 

Limited by the signal-to-noise ratio (S/N) for the 
individual-night spectra, there are degeneracies between sev¬ 
eral pairs of spectral components, including: (a) the AGN 
power-law continuum and the host galaxy; (b) Fe II emission 
and the broad He II line; (c) Fe II emission and the coronal 
lines. Thus, we first fit the high-S/N mean spectrum with all 
parameters set free. We then fit each individual-night spec¬ 
trum with the values of some parameters fixed to those ob¬ 
tained from the fitting of the mean spectrum. The details of 
each spectral component, and the fitting parameters, are de¬ 
scribed in the following subsections. 

9 http://ned.ipac.caltech.edu/ 

10 Note that the pseudo-continuum here is not defined in the ordinary way, 
in which the host galaxy component is not included. 


3.1.1. AGN Power-law Continuum 

A single power law is used to describe the featureless AGN 
continuum. It has two parameters: the flux density at 5100 
A (Fagn)j an d the spectral index (a, defined as j\ oc A“). 
There is some degeneracy between the power-law continuum 
and the host galaxy. A larger a or a higher-flux galaxy (,Fg a i) 
both make the total spectrum redder (note that the unphysi¬ 
cal change in the color of the observed spectrum caused by 
weather or differential atmospheric refraction is avoided in 
our observation, as our flux calibration by the comparison star 
provides differential spectroph otom etry in each wavelength 
bin of the spectrum; see Section [T2l i. 

For our observations, we find that the relative flux of the 
galaxy component is the main reason for the change in the 
total-spectrum slope. Figure |T|illustrates this point. The right 
panel shows the spectrum of Mrk 382 taken at JD 2456299, 
just one day after the date of the spectrum shown in the left 
panel. Apparently, the flux at 5100 A is ~30% higher than 
the day before, and the color is redder. However, the differ¬ 
ence between the V -band magnitudes of the two nights is only 
~0.01 mag. We have therefore adopted an approach based on 
the assumption that the absolute flux of the host must be con¬ 
stant but the relative flux inside the slit can vary depending 
on the observing conditions. This is illustrated in detail in 
Appendix [A] The conclusion about the AGN continuum is in 
line with several, but not all studies discussing the relation 
between a and luminosity in AGNs. Most recently, IZhand 
(12013ft investigated this issu e using the reverbe ration mapping 
data of the 17 quasars from iKaspi et al.l d2000l) . They found 
no strong dependence of a on the variability of the luminosity. 
Thus, in the fitting of each individual-night spectrum, we fix 
the value of a to agree with the best fit to the mean spectrum, 
and let Fagn and F ga i free. Comparing the left and right pan¬ 
els of Figure |T] only the galaxy component changes. In this 
way, the resultant Fagn light curve matches the V -hand light 
curve well (see Appendix[A]i. Our fitting is consistent with no 
changes in spectral index as a function of luminosity. 

3.1.2. Few Emission 

There are several optical Fe II templates available in the lit^ 
erature . A mong them, two templates , from lBoroson & Greenl 
d 1992b and IVeron-Cettv et al.l (12004ft . are most widely used, 
both constructed from the spectrum of the NLS1 galaxy 

I Zw 1 . IBarth et al.1 (12013ft find that the template from 
IVeron-Cettv et al.l (l2004li ~ yields inconsistent He I emission 
lines; significant broad He I A4922 and A5016 emission lines 
are needed while He I A4471 is constrained to have zero flux 
by the fitting. We compared these two templates and found 
that the one from lBoroson & Greenl (11992ft gives better fitting 
as judged by the smaller reduced x 2 - No broad He I A4922 
and A5016 emission lines are needed, which is consistent with 
the zero flux of He I A4471. Introducing the two broad He I 
lines will add more parameters to the fit and will introduce ad¬ 
ditional degeneracy between He I A4922 and H/3. Given th is, 
we chose the Fe II template from lBoroson & Greenl d 1992ft in 
this paper. 

The Fe II template is convolved with a Gaussian function 
to be scaled, broadened, and shifted (see IHu et alj|2008l for 
details). Three parameters, the flux (Fp e ), the full width at 
half-maximum (FWHMp e ), and shift (Vp e ), are used for Fe 

II emission in the fitting. We let all the three parameters to 
be free in the fitting of each individual-night spectrum. Note 
that, Fp e is defined as the flux of the integrated Fe II emis- 
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Figure 1. Individual-night spectra and multi-component model fitting for Mrk 382 taken at JDs 2456298 (left) and 2456299 (right). For each night, the top 
panel shows the Galactic extinction and redshift corrected spectrum (green or black for the parts in/out of the fitting window), and best-fit model (red). The model 
includes a pseudo-continuum (blue, composed by AGN power-law continuum, Fe II emission, and host galaxy) and several emission lines (FI/3 in magenta, broad 
He II A4686 in cyan, and narrow emission lines in orange). The bottom panel shows the residuals. Note the large change in the host galaxy intensities for the two 
successive nights, which is not real but caused by our method of observation and calibration (see Appendix lAlfor details). 


sion betw een 4434 and 4684 A fro m the best-fit Fe II model, 
following lBoroson & Greenl (119921) . 

3.1.3. The H/3 Line 

The spectral resolution of our spectra is rather low (~500 
km s _1 ), and in addition our objects are selected to have nar¬ 
row H/3 emission lines. Thus, it is hard to decompose the H/3 
emission line and remove the narrow component that comes 
from the narrow-line region. Considering that the contribu¬ 
tion of the narrow component is weak and supposed to be 
constant over reverberation timescales dPeterson et al.l[2(Jl3 v ). 
we treat the entire H/3 emission line as one component. The 
entire H/3 profile is modeled by a Gauss-Hermite function 
(Ivan der Marel & Franxl fl~993b . All the five parameters of 
the Gauss-Hermite function are set free in the fitting of each 
individual-night spectrum. The flux of the line is calcu¬ 
lated from the best-fit model. 

3.1.4. Narrow Emission Lines 

Besides the strong [O III] AA4959, 5007, there are many 
other narrow emission lines in the wa velength range of our 
fit spectrum (Ivanden Berk et all 120011) . We identify narrow 
He II A4686, He I A4471, and several high-ionization forbid¬ 
den coronal lines. In some objects, these narrow lines are so 
strong that their contamination is non-negligible to the Fe II 
measurement. Figure [2] shows an individual-night spectrum 
of Mrk 335 and our fit. The strong coronal lines include in 
this case [Fe vn] A5158, [Fe vi] A5176, [N i] A5199, [Ca v] 
A5309, [Fe VII] A5721, and [Fe VII] A6086. Adding these 
lines significantly improve the fit. 

The intensities of the narrow lines differ from one object 
to the next. For each object, we first identify narrow lines in 
the mean spectrum by testing the goodness of fit. Only lines 
that are identified in the mean spectrum are included in the 
fitting of individual-night spectra. Each narrow emission line 
is modeled by a Gaussian. These lines are coming from the 
narrow-line region and hence do not vary on the campaign 
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Figure 2. An example of individual-night spectrum and model fitting for 
Mrk 335. The spectrum, model, and residuals are plotted in the same manner 
as that in Figure^ Note the strong narrow He II and coronal lines labeled by 
ion and wavelength. 


time scale. Thus, we constrain them to have the same velocity 
width and shift with those of the [O in] A5007 line. The flux 
ratio relative to [O ill] A5007 is also kept constant, as given 
by the best fit to the mean spectrum (the flux ratio of [O in] 
A4959 to A5007 is fixed to the theoretical value of 1/3), in 
order to avoid the degeneracy between these lines and Fe II. 
Only the flux, velocity width, and shift of [O III] A5007 are 
set free. 


3.1.5. Broad He II Emission Line 
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The broad He II A4686 emission line is strong in the spec¬ 
trum of some objects in our sample. Figure[2]shows an exam¬ 
ple. The line, plotted in cyan, is much broader than H/3, and 
strongly blueshifted. Such a line profile is common for high- 
ionization lines in NLS1 s and consistent with that of the ultra¬ 
violet (UV) He II A1640 line in SDSS spectra ([Richards et al.l 
1201 ll) . It is even more prominent in the root-mean-square 
(rms) spectra of our objects, indicating large variability sim¬ 
ilar to the r esults in previous He II reverberation mapping 
studie s (e.g., iBentz et aI1l2010t iBarth et aTll20111 iGrier etaH 

I20l2h . 

The broad He II line is heavily blended with Fe II emis¬ 
sion. The high-ionization lines in AGN are often asymmetric 
([Richards et al.l l201lb . but there is no other high-ionization 
line in the wavelength range of our spectra to constrain the 
profile of He II line. We first fit the line profile in the mean 
spectrum using a single Gaussian. We then kept the width 
and the shift in the individual-night spectra but leave the line 
intensity as a free parameter. 

The intensity variation of the He II lines in our sample will 
be discussed in a forthcoming publication. Here we only note 
that the procedure described above seems to be consistent 
with the observations, and, in general, the lag of this line is 
much shorter than the H/3 lag. 

3.1.6. Host Galaxy 

Following B arth et afl (120131). we use sing l e simp le stel¬ 
lar population models from lBruzual & Chariot! (120031) as tem¬ 
plates for the galaxy component. For most objects in our sam¬ 
ple, the instantaneous-burst model with an age of 11 Gyr and 
solar metallicity (Z = 0.02) provided sufficiently good fit to 
the mean spectrum and consistent flux ratio compared with 
the measured Hubble Space Telescope (HST ) images (see Pa¬ 
pers I and II). In a few cases (Mrk 335, Mrk 142, and Mrk 42), 
this template gives a flux that is much larger than that derived 
from the HST image. The template with 11 Gyr and Z = 0.05 
provides better results and was adopted in these cases. 

Like the Fe II template, the galaxy template is also con¬ 
volved with a Gaussian to be scaled, broadened, and shifted. 
As described previously, and detailed in Appendix [A] the rel¬ 
ative flux /\, a | is free to vary in the fitting of individual-night 
spectra. 

3.2. Fitting Results 

Table |T]lists the values of the parameters fixed in the fitting 
of individual-night spectra for all the objects in our sample. 
Columns (5) to (12) list the narrow emission lines included 
in the fits. As explained, their velocity widths and shifts are 
constrained to be the same with those of the [O III] A5007 
line while the relative intensity ratios are kept as measured 
from the mean spectrum. The last column gives the galaxy 
template used. In the table, a blank entry means that spectral 
component is not included for the object. 

Figure [3] shows examples of fittings to individual-night 
spectra for the eight objects not presented so far. The nota¬ 
tions and colors are as the same as those in Figure|T]for Mrk 
382. Notes on two objects with unusual treatments are as fol¬ 
lows. 

IRAS 04416+1215. This object has the highest redshift 
(0.089) in our sample. From the images taken for slit cen¬ 
tering, the point spread function of the object is as broad as 
that of the comparison star, meaning that there is a negligi¬ 
ble galaxy contribution. For this object, the comparison star 


is fainter than the AGN, especially in the red part. We found 
that the shape of the spectrum of this object is not as well cali¬ 
brated by the comparison star as in the other objects. Thus, we 
fit in a relatively narrow wavelength window (4430-5550 A), 
but covering the two bumps of Fe II emission on either side of 
H/3. We also let the spectral index a of the AGN power-law 
continuum free to vary in individual-night spectra to compen¬ 
sate for the difference in spectral shape caused by the calibra¬ 
tion. As shown in Figure^ the spectrum can be fitted well by 
including only a power law, Fe II emission, H/3, [O III], and a 
few narrow lines. 

IRAS F12397+3333. As discussed in Paper I, the spec¬ 
trum of this AGN is probably affected by host galaxy extinc¬ 
tion. The observed spectra (see Figure 4 of Paper I for an 
example) have a very red color, and are badly fitted by our 
spectral model. We perform an extinction correction after de- 
redshifting, assuming the Galactic extinction law and Ay = 
1.44 mag estimated from the Balmer decrement (see details 
in Paper I). This dereddened spectrum is fit well with a typi¬ 
cal spectral index as shown in Figure [3] 

3.3. Line-profile Measurements 

For each individual-night spectrum, the fluxes, FWHMs, 
and velocity shifts (with respect to [O III]) of Fe II and H/3 
are calculated from the best-fit model. The means of these 
properties are considered to be the measurements of the line 
profile, except for the width of H/3 (FWHMh/ 3 ; see below). 
The standard deviations on those means are used as the un¬ 
certainties. 

FWHMh /3 is underestimated in the fitting, as we do not in¬ 
clude the narrow component of H/3 in our model (see Section 
Ed- Thus, we use a different method to estimate this width. 
We add a Gaussian to our model to represent the narrow H/3 
component. The velocity width and shift are constrained to 
be the same as those of [O III], We then fit the mean spectra 
twice, once assuming the flux of the narrow H/3 line is 10% 
of the flux of [O III] A5007 and once assuming a flux ratio of 
0.2. The width of the broad H/3 obtained with flux ratio set to 
0.1 is adopted as the FWHM^. The uncertainty is obtained 
from the fit with a flux ratio set to 0.2 and the original fit that 
assumed no narrow H/3. This method essentially resembles 
the one used in Papers I and II, while the spectral model here 
is more sophisticated. 

Table[2]lists the measurements of the 10 objects. The ins tru- 
mental broadening (FWHM ~ 500 km s -1 . Section ^. 1 .31 ) has 
been taken into account in the listed FWHMs. Note that the 
H/3 of IRAS 04416+1215 has large velocity shift with respect 
to [O III] A5007, and its profile shows no feature of a narrow 
component with the same shift of [O III], Adding a narrow 
component makes the width of the broad H/3 even narrower. 
Thus, the FWHM H/9 of IRAS 04416+1215, as other proper¬ 
ties, was obtained from the measurements of individual-night 
spectra. 

3.4. Light-curxe Measurements 

Our fitting successfully reduces the scatter in the light 
curves due to the influence of the host galaxy contamination; 
thus, Fagn represents the real AGN continuum much better 
than the simply integrated 5100 A flux OF5100X as shown in 
Appendix lAl The light curves of the emission lines are also 
generated directly from the best-fit values of the correspond¬ 
ing parameters obtained from the fits of individual-night spec¬ 
tra. The errors of the fluxes given by the fitting are not large 
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Table 1 

Parameters Fixed in the Fitting of Individual-night Spectrum 


Object 

(1) 

Power Law 

a 

<Jx °C A“) 
( 2 ) 

Broad He II A4686 



Narrow Emission Lines 



Galaxy model 

(13) 

FWHM ' 1 

(km s' 1 ) 

(3) 

Shitt 

(km s' 1 ) 

(4) 

He I 

A4471 

(5) 

He II 

A4686 

(6) 

|Fe VIIJ 
A5158 
(V) 

|Fe VIJ 
A5176 
( 8 ) 

IN I] 

A5199 

(9) 

LCa V] 

A5309 

( 10 ) 

|Fe VIIJ 
A5721 
( 11 ) 

[Fe VIIJ 
A6086 
( 12 ) 

Mrk 335 

-1.45 

7082 

-620 

0.042 

0.115 

0.043 

0.032 


0.028 

0.057 

0.079 

1 lGyr z05 

Mrk 1044 

- 2.01 

5161 

-942 

0.216 


0.117 

0.219 



0.205 

0.225 

1 lGyr_z 02 

IRAS 04416+1215 

free b 






0.095 

0.042 





Mrk 382 

-1.96 

5797 

-1322 

0.028 

0.074 




0.021 

0.033 

0.021 

llGyr z 02 

Mrk 142 

- 2.11 

4824 

-757 

0.041 

0.131 

0.055 

0.116 



0.082 

0.095 

1 lGyr z05 

MCG +06-26-012 

-0.84 

5191 

-927 

0.063 

0.096 





0.046 

0.061 

1 lGyr z 02 

IRAS F12397+3333 c 

-2.28 

5506 

-334 


0.034 

0.013 


0.018 


0.008 

0.010 

1 lGyr z 02 

Mrk 42 

-0.70 

1472 

-55 

0.069 

0.093 







1 lGyr z05 

Mrk 486 

-0.76 

4877 

-724 


0.171 

0.027 



0.067 

0.050 

0.140 

llGyr z02 

Mrk 493 

-0.91 

4054 

-1580 

0.132 

0.261 







llGyr z02 


Note. — Fixed parameters in the fitting of individual-night spectra. Columns (2) to (4) list the absolute values of the parameters. Co lumns (5) to (12) list 
the relative intensity ratios with respect to [O III] A5007 for the narrow emission lines. Column (13) lists the galaxy template from tBruzual & Chariot! (20031) . 
• • • means that the component is not included in the fitting for the specific object. 


a The FWHMs listed are after instrumental broadening correction. 

b The spectral index of the power-law continuum is free to vary in the fitting of the IRAS 04416+1215 spectra, see text for details. 
c Intrinsic host galaxy extinction is assumed in the fitting of IRAS F12397+3333. See text for details. 


Table 2 

Line-profile Measurements 


Object 


Fe II 



H/3 



Mux 

JbWHM 

Shitt 

Flux 

FWHM- 

Shitt 


( 10' 15 erg s ' 1 

cm -2 ) (kms -1 ) 

(km s' 1 ) 

(10 15 erg s ' 1 cm' 2 ) 

(km s' 1 ) 

(km s' 1 ) 


Mrk 335 

253 

± 

9 

1947 

± 

143 

49 

± 

37 

661 

± 

22 

2096 

± 

170 

-1 

± 

16 

Mrk 1044 

369 

± 

11 

866 

± 

25 

76 

± 

10 

380 

± 

16 

1178 

± 

22 

12 

± 

12 

IRAS 04416+1215 

292 

± 

9 

1313 

± 

50 

496 

± 

21 

149 

± 

4 

1522 

± 

44 b 

241 

± 

29 

Mrk 382 

27 

± 

4 

1326 

± 

234 

13 

± 

115 

39 

± 

2 

1462 

± 

296 

-44 

± 

35 

Mrk 142 

87 

± 

5 

1512 

± 

69 

-25 

± 

37 

78 

± 

5 

1588 

± 

58 

-101 

± 

28 

MCG +06-26-012 

41 

± 

4 

1155 

± 

70 

-21 

± 

41 

41 

± 

4 

1334 

± 

80 

-29 

± 

23 

IRAS F12397+3333 

550 

± 

28 

1748 

± 

78 

36 

± 

40 

405 

± 

18 

1802 

± 

560 

-50 

± 

26 

Mrk 42 

67 

± 

3 

787 

± 

16 

105 

± 

19 

57 

± 

2 

802 

± 

18 

87 

± 

19 

Mrk 486 

186 

± 

5 

1790 

± 

88 

95 

± 

33 

346 

± 

12 

1942 

± 

67 

-46 

± 

9 

Mrk 493 

102 

± 

3 

780 

± 

9 

171 

± 

9 

92 

± 

3 

778 

± 

12 

126 

± 

13 


Note. — Fluxes, FWHMs, and velocity shifts of Fe II and H ft. Except for the HO FWHMs, the listed values are the means and 
standard deviations obtained from the measurements of individual-night spectra. All listed FWHMs include a correction due to 
instrumental broadening. The velocity shifts are with respect to [O III], 


- H O FWHMs are estimated from the mean spectrum of each object taking into account the narrow H O component, as explained in 
the text. 

b H/3 FWHM for IRAS 04416+1215 is obtained from individual-night spectra. See text for details. 


enough to account for the scatter in the fluxes of successive 
nights, an additional systematic error is estimated for each 
light curve (as in Paper I). This systematic error is added in 
quadrature to the fitting error for the calculations of variability 
amplitudes and time lagsbelow. Note that our treatment is dif¬ 
ferent from that of lBarth et all (120131 ). who measured the light 
curve of H/3 by integrating the continuum-subtracted spectra. 

4. LIGHT-CURVE ANALYSIS AND RESULTS 

The left columns of Figure 0] and Figure 0 show the light 
curves of Fagn, Fh/s, and Fp e for nine objects with reliable lag 
measurements, and the remaining one (Mrk 42), respectively. 
In this section, we will calculate the variability amplitudes and 
the reverberation lags for the H/3 and Fe II emission lines. We 
then compare the H/3 time lags with those presented in Papers 
I and II. 

4.1. Variability Amplitudes 

We use the quantity F val defined in lRodrfguez-Pascual et all 
ill9971) to represent the variability amplitude. This quantity is 


an estimate of the intrinsic variability over the errors. The un¬ 
certainties are calculated following lEdelson et al.ld2002l) . The 
results for Fe II and H/3 are listed in Table[3] For the variabil¬ 
ity amplitude ratio of Fe II to H/3, the range is from about 
0.6 (Mrk 486) to 1.2 (Mrk 42), except for Mrk 382 which has 
the largest host galaxy contamination. On average (except 
for Mrk 382), the Fe II value is about 1 0% sm aller, consis¬ 
tent with previous results (e.g. JVestergaard & Peterso~nll2005i : 
iBarth et alMIl) . 

4.2. Reverberation Lags 

The time lags between the AGN continuum variations 
(Tagn) and the emission lines (Fp e and Fhh) are measured 
from the cross-correlation functions (CCFs) for the rele¬ 
vant light curve s. We use the interpolat i on cross-correlation 
funct i on (ICCF: [G askell & Sparkel[l986l : lGaskell & Petersonl 

1987; Whi te & Petersonll 1 9941) method to calculate the CCF, 
and adopt the centroid of the CCF, above 80% of the 
peak value (r max ) as the time lag dKoratkar & Gaskelllll991l : 
iPeterson et alJ l2004). The uncertainty in the time lag mea- 
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Figure 3. Examples of individual-night spectra and model fittings for the other eight objects. For each object, the spectrum, model, and residuals are plotted in 
the same manner as those in Figure [Tj 


Table 3 

Variability Amplitudes 


Object 

Fvar 


Fen Fl/3 


Mrk 335 

3.1 

± 

0.3 

3.0 

± 

0.3 

Mrk 1044 

2.6 

± 

0.3 

3.7 

± 

0.4 

IRAS 04416+1215 

2.1 

± 

0.3 

2.0 

± 

0.3 

Mrk 382 

11.5 

± 

1.2 

4.1 

± 

0.4 

Mrk 142 

5.5 

± 

0.5 

6.6 

± 

0.5 

MCG +06-26-012 

8.1 

± 

1.2 

9.2 

± 

1.2 

IRAS F12397+3333 

4.4 

± 

0.6 

4.1 

± 

0.5 

Mrk 42 

3.6 

± 

0.7 

2.9 

± 

0.4 

Mrk 486 

2.0 

± 

0.5 

3.4 

± 

0.4 

Mrk 493 

2.1 

± 

0.7 

3.1 

± 

0.5 


Note. — The listed values are in percentage. 


surement is estimated from the cross-correlation centroid dis¬ 
tribution (CCCD) given by random subset selection/flux ran- 
domization Monte C arlo realizations (Mao z & Netzerl fT989t 
iPeterson et al.ll 1998b . 

All objects in our sample, except for Mrk 42, have reli¬ 
able time lag measurements for both Fe II and H/3. The right 
columns of Figure [4] show the results of the CCF analysis. 
For each object, the top-right panel shows the autocorrelation 
function (ACF) of the Fagn light curve, which is shown in the 
top-left panel. The two lower panels in the right column show 
the CCFs (in black) for H/3 (middle-right) and Fe II (bottom- 
right) with respect to Fagn light curves. The blue histograms 
are the corresponding CCCDs. Most CCCDs have a rather 
symmetric profile, except those for IRAS 04416+1215. 

Table |4]lists the time lags of Fe II (rp e ) and H/3 (th/j), their 
uncertainties, and the corresponding r max , for the nine objects 
with reliable lag measurements. The listed time lags and un- 
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Figure[3] Continued. 


Rest Wavelength (A) 


certainties are in the rest frame after time-dilation correction. 
Both time lags for IRAS 04416+1215 have uncertainties of 
highly unequal upper and lower limits, as the sequence of 
their asymmetric CCCDs. For other objects, the time lags 
are well determined by the single-peak CCFs with high r max 
and symmetric, narrow CCCDs. 

The detection rate of Fe II time lag in our sample is ex¬ 
tremely high (9 in 10), compared with previous reverberation¬ 
mapping experiments. This is mainly attributed to the com¬ 
mon feature of strong Fe II emission in the SEAMBH sample, 
the high-cadence observation, and sophisticated spectral fit¬ 
ting. The failure to obtain reliable 7p e and Tup in Mrk 42 is 
caused by the Fagn light curve, which has large scatter and 
no clear structure (see Figure 0). However, the light curves of 
Fp e and Fh /3 show hints for a similar structure which, unfortu¬ 
nately, is not enough to establish a time lag. 

The clear reverberation of the Fe II emission in response to 
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Table 4 

Time Lag Measurements 


Object 

Fe 11 


H/3 

r max 

TTe 

r max 

7H/3 

Mrk 335 

0.48 


0.70 

8.7+i'g 

Mrk 1044 

IRAS 04416+1215 

0.48 

0.61 

13.9^ 

12 - 6!$ 7 

0.62 

0.60 

I0 -5j7 

17} 72+0.9 

Mrk 382 

0.40 

23-8-60 

0.47 

7-5+11 

Mrk 142 

0.78 

7 6 +1 - 4 

0.84 

7.9$ 

MCG +06-26-012 

0.86 

22 4 +9,3 
z ' z, ‘ -6.3 

0.92 

24.0^-g 

IRAS F12397+3333 

0.63 

10 . 6 +™ 

0.71 

9.7+jj 

Mrk 486 

0.77 

17.3TI 

0.79 

23.7$ 

Mrk 493 

0.63 

11 9 +36 

1 1 **-6.5 

0.81 

11-6 $ 


Note. — Time lags are in the rest frame, and in units of days. 


the continuum supports the hypothesis that the Fe II emission 
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Figure 4. Light curves and cross-correlation functions for the nine objects with reliable time lag measurements. For each object, the left column show the light 
curves of the AGN (top-left), H (3 (middle-left), and Fe II (bottom-left). The error bar of each flux point is only the fitting error. For each light curve, an additional 
systematic error is estimated from the scatter in the fluxes of successive nights, and plotted as the error bar with terminals in the comer of the panel. The right 
column show the autocorrelation function of the AGN light curve (top-right), cross-correlation functions for H/3 (middle-right) and Fe II (bottom-right) with 
respect to the AGN continuum. The blue histograms in the middle-right and bottom-right panels are the corresponding cross-correlation centroid distributions. 


originates from photoionized gas in the BLR. The high detec¬ 
tion rate indicates that this is prevalent in NLS1, which are 
AGNs with high Lboi/^Edd- It is possible that the origin of Fe 
II emission in AGNs with low Lboi/TEdd is different, but there 
is little evidence to support this claim. In fact, iBarth et ail 
(120131) has already detected the reverberation of Fe II in two 
broad-line Seyfert 1 galaxies with H/3 widths of FWHM > 
4000 km s -1 and much lower Lboi/TEdd- Comparing the time 
lags of Fe II and H/3, along with their velocity widths, pro¬ 
vides some hints about the location and geometries of the two 
emission regions. We will discuss these ideas in Section 15721 
below. 

4.3. Hf3 Lag Comparison 

In Paper I and II, the reported time lags of H/3 are obtained 
from the light curves measured by the more traditional inte¬ 
gration method, without taking into account the contamina¬ 
tion by the host galaxy and the narrow emission lines. Figure 


[ 6 ] shows a comparison between the previously reported Tup 
based on the integration method and the new results given by 
the fitting method reported here. The dashed diagonal line 
denotes the 1:1 ratio and the region between the two dotted 
lines gives the 0.1 dex deviation from the line. For seven 
out of nine objects, the differences between Tup given by the 
two methods are less than 0.1 dex. The two exceptions are 
IRAS 04416+1215 and Mrk 1044. In Paper II we show that 
using the integration method, we cannot obtain a significant 
th /3 for IRAS 04416+1215. The lag for Mrk 1044 reported 
in that paper is about a factor 2 shorter than obtained here, 
with a very large uncertainty (4.87] 7 ). This value is consis¬ 
tent with the new one presented here within the uncertainties. 
In addition, the uncertainty on Tup in Mrk 493 given by the 
integration method is very large with a negative lower limit 
(12.27]g 7 ). The fitting method yields much smaller uncer¬ 
tainty with basically the same Tup. This is a direct result of 
the much smoother H/3 light curve of Mrk 493 given by the 
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Figure[4] Continued. 




Figure[4] Continued. 


fitting method (see Figure [4j and also Figure 2 of Paper II). 
Thus, the time lags of H/3 given by the two methods are con¬ 
sistent with each other and r H /-j is better defined by the fitting 
procedure in some of the cases. 

5. DISCUSSION 

5.1. The Radius-Luminosity Relationship for Fe II 

A simple theoretical expectation based on photoioniza¬ 
tion is that /?blr oc L a with a ~ 0.5. The Rblr-L rela¬ 
tion for Hij has been compare d with such predictions in sev¬ 
eral earlier publications (e.g., Kaspi el al. 2000: B entz et all 
20091 12013L and refere nces therein). Regarding Fe II lines, 
Chelouche et ali (120141) provide a tentative R^lr-L relation 
for the first time, from a small inhomogeneous sample of six 
AGNs. Among them, two objects are reported bv lBarth et all 
(120131) usin g the standard spect ral fit ting method, four others 
come from Ra fter et aD (120131) and IChelouche et all (120141) 
using the MCF scheme of IChelouche & Zuckerl(l20T3 ). Here, 
we revisit this issue by more than doubling the size of the 
sample by adding our nine newly_measured sources. For a 
proper comparison with the lChelouche et aU(l2014l) results we 
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Figure 5. Light curves of the AGN (top panel), H/3 (middle), and Fe II 
(bottom) for Mrk 42. Note the similarity of the light curves of H/3 and Fe II. 



running) (Days) 

Figure 6 . A comparison of the H/3 time lags obtained by simple integration 
(Papers I and II) with those obtained by fitting (this work). The differences 
are typically less than 0.1 dex (dotted lines) except for IRAS 04416+1215 
(not plotted because its th /3 cannot be defined by integration) and Mrk 1044 
(very large uncertainty of the integration method). Mrk 493 also has very 
large uncertainty when measured by the integration method but the measured 
th /3 are very similar. 


used our best observed fluxes and uncertainty and a cosmol¬ 
ogy with //(> = 70 km s -1 Mpc~', D m = 0.3, and 17 a = 0.7. 

Figure [7] shows TF e vs. the 5100 A luminosity for the nine 
objects in our sample (black dots) alongside the earlier re¬ 
sults (kindly provided by D. Chelouche). The green squares 
are those obtained by the MCF scheme (including three with 



Figure 7. Fe II Rblr-L relation. The nine objects in our sample are plotted 
as black dots, the three objects with their names labeled have relatively large 
differences between and rp e - Measurements from previous studies are 
plotted as blue triangles and green squares (see text for details). The solid line 
is the fit to the entire data set. The dotted and d ashed lines are the Rblr-L 
relations in I Chelouche et all <20141) for Fe II and IBentz et all <20091) for H/3, 
respectively. 

insignificant results, see Figure 4 of IChelouche et all 120141 
for details). The blue triangles are the two objects from 
iBarth et akl (120131) for which we used the centroid of the pub¬ 
lished CCFs for Fe II and the F-band light curves. Using the 
FITEXY method dPress et al.ll 19921) to fit all the 18 objects we 
find: 

log(/?Fe i. blr) = (—22.0± l.l)+(0.54±0.02)log(AL A (5100 A)), 

0 ) 

which is plotted in a solid line . This line deviates cons iderably 
from the relation presented in lChelou che et all (120 14l) (shown 
in a dotted line), but close to that of lBentz et al .1(1200 9) for H/3 
(dashed line). This result is expected since most of the ob¬ 
jects in our sample show roughly equal Tp e and tu /3 (see Table 
EU, except for three objects (Mrk 335, Mrk382, and Mrk486), 
while the previous studies have rp e greater than t H/ 5 . Our sam¬ 
ple by itself shows no correlation at all between rp e and the 
5100 A luminosity. This is not surprising given the relatively 
small range in luminosity considered by us and the fact that 
for this sample, the correlation between tha and the contin¬ 
uum luminosity is also very weak. We suggest that this lack 
of correlation is related also to the large Lboi/^Edd of the ob¬ 
jects in our sample, a topic which we discuss in great detail in 
Paper IV (Du et al., submitted). 

5.2. Comparison between Fe II andF[f3 

As shown in Table [4] rp e is roughly equal to r ( i; j (count¬ 
ing the uncertainties in both lags) in six of our objects, longer 
in two (Mrk 335 and Mrk 382) and somewhat shorter in one 
(Mrk 486). A common feature of the three objects with un¬ 
equal lags is the relatively low Fe ll/H/3 intensity ratio. Fig¬ 
ure [8] shows the plot of rp e /rH /3 vs. FpJFnp (defined as Rp e ). 
The three objects with different time lags are labeled by their 
names and locate in the region of Rp e < 1. The six other ob¬ 
jects show Rp e > 1. We also mark the approximate positions 
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Figure 8. Time lag ratios of Fe II to H (3 vs. intensity ratios of Fe II to H/3. 
Objects in our sample are shown as black dots. The three objects labeled by 
names have different time lags of Fe II and H/3, and all have low intensity 
ratio of Fe II to H/3 . Bl ue tria ngles show the approximate positions of the two 
objects in IBarth et al.l (2013D (intensity ratios are estimated by eye). The two 
points are slightly shifted horizontally for clarity (See text for details.) 

of the two objects measured bv lBarth et al.l (120131)1^1 Based 
on eye estimates of the published light curves shown in their 
paper, both objects have Rp e ~ 0.75. Their positions in Figure 
[5] are marked around 0.75 (and horizontally shifted for clar¬ 
ity). They both have different time lags than H/3 and Rp e < 1, 
consistent with the results of our sample. 

We note, again, that our sample is highly biased toward 
AGNs with high Lboi/^Edd and thus high Rp e . In an unbiased 
AGN sample, the fraction of sources with Rp e < 1 is much 
larg er (see, e.g.. Figu re 7 of ISu ent ic et al.l l2000 and Figure 

I of lShen & Holl2014t) . iHu et al. ( 120081 ) analysed the proper¬ 
ties of a large number of low-redshift AGNs. They found that 
the velocity width of Fe II is systematically narrower (^3/4) 
than that of H/3, which may suggest, given Keplerian veloci¬ 
ties, that in those sources tv c is longer than r H . 3 . Our present 
sample contains only 10 sources which is too few to test, sys¬ 
tematically, any of these ideas. 

6 . SUMMARY 

We provide new Fe II measurements for 10 NLSls and re¬ 
port statistically significant time lag measurements for nine 
of these sources. All the observed NLSls are suspected to 
be SEAMBHs. This more than doubles the number of AGN 
showing measurable Fe II time lags. Our time lag measure¬ 
ments are based on high-cadence, high-S/N measurements at 
the Yunnan observatory, and on a new, sophisticated, fitting 
analysis that takes into account the uncertainties caused by 
the apparent flux change of the host galaxy and several nar¬ 
row emission lines. We demonstrate, by a careful comparison 
with earlier measurements of th/ 3 , that the method can con¬ 
siderably improve the accuracy of the time lag measurement. 

11 The Fe 11 flux in IBarth et all ( 120131) is defined as the integrated flux of Fe 

II between 4400 and 4900 A, while our definit ion refers to the range bet ween 
4434 and 4684 A. For the Fe II template of IBoroson & Greenl 119921) . Fp e 
defined here is about 3/4 that in Barth et al. 1 2QT3). 



The main findings reported in this paper can be summarized 
as follows: 

1. All 10 objects presented here show Fe II variations with 
an amplitude of a few to ten percent. On average, this 
is about 10% smaller than the variability of the H/3 line. 

2. Reliable Fe II reverberation time lags with respect to 
the AGN continuum are detected in nine objects, con¬ 
firming the suggestion that the Fe II emission originates 
from photoionized gas. 

3. Combining the new reverberation mapping results 
with those in previous work, shows a clear radius- 
luminosity relationship for Fe II that is similar to the 
one for H/3. However, our sample by itself shows no 
such correlation due to the large intrinsic scatter over a 
small luminosity range. 

4. The difference in the time lags of Fe II and H/3 depends 
on the intensity ratio of Fe II to H/3 (Rp e ). The time 
lag of Fe II is roughly equal to that of H/3 in all the six 
objects with Rp e > 1. The Fe II time lag is longer in 
Mrk 335 and Mrk 382, those objects with Rp e < 1, and 
shorter in Mrk 486 with Rp e < 1. 
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APPENDIX 

A. HOST GALAXY FLUX CALIBRATION 

Our flux calibration method provides the relative flux of the 
object with respect to that of the local comparison star. For 
two point sources kept in a line parallel to the slit, the fractions 
of light loss due to seeing, differential atmospheric refraction, 
and mis-centering are identical. Thus this method gives ac¬ 
curate relative flux calibration for the spectral components of 
the AGN, including the featureless power law (//\g\) and the 
broad emission lines. This is not the case for the extended 
host galaxies whose flux relative to the comparison star inside 
the slit can vary much more due to seeing variations and mis- 
centering. In addition, different parts of the galaxy may be 
observed each time. As a result, the derived flux of the host 
galaxy is F g!l \ = /caiFgatabs, where F ga i ja b s is the absolute flux 
(which is constant), and / ca 1 is a factor accounting for all the 
effects of varying observing conditions in the flux calibration 
procedure. / ca 1 could change from one exposure to the next, 
introducing an additional uncertainty into the integrated 5100 
A flux 7*5100 = F A GN + /calFgal,abs- 

In Papers I and II, we estimated the flux of the host galaxies 
in the spectral extraction aperture using archival HST images 
for eight objects in our sample. The resultant relative fluxes, 
(/•„ai///\GN). range from ~0.2 in Mrk 335 to > 1 in Mrk 382. 
Thus, the apparent change of flux of the host galaxy is not 
negligible for those objects with strong host galaxy contribu¬ 
tion. For Mrk 382, the measurement of the AGN continuum 
was badly affected by this uncertainty, which forced us to use 
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Figure 9. Light curves of the integrated 5100 A flux (F 5100 , top panel), fit¬ 
ted host galaxy (/cai^gai^bs described in Appendix [A] upper-middle panel), 
fitted AGN power-law continuum (Fagn> lower-middle panel), and V-band 
magnitude (bottom panel, see also Paper II). Note that F 5100 is heavily con¬ 
taminated by the host whose apparent flux is varying by f ca \, and the Fagn 
recovers, accurately, the intrinsic AGN variability as judged by a comparison 
with the V-band observations. 


the V-band photometry instead. The new procedure adopted 
here (see the main text) enables us to solve for / ca i-Fgai,abs for 
each individual-night spectrum. The detailed fitting of the 
various spectral components considerably reduces the noise 
in the various line and continuum light curves. In particular, 
for Mrk 382, we can now recover the 5100 A AGN continuum 
variations to a much better accuracy, which is evident by the 
good agreement with the V-band photometry. 

The big improvements due to the spectral fitting procedure 
are shown in Figure 0 which provides more information on 
the process for Mrk 382. The top panel shows the light curve 
of F 5100 , which is measured as in Papers I and II, except for the 
inclusion of Galactic extinction correction and de-redshifting. 
The scatter in this light curve is very large, and hence it was 
replaced in Paper II by the V-band light curve. The upper- 
middle and lower-middle panels show the results of the spec¬ 
tral fitting for F ga i and Fagn> respectively. The light curve 
of Fagn is basically identical to that of the V-band (bottom 
panel). Clearly, the unrealistic high F 5100 around JD 2456300 
was caused by the large contribution of the host galaxy caused 
by seeing and guiding fluctuations. On other JDs, the host 
contamination is relatively small, but still strong enough to af¬ 
fect the variability of Fagn- e.g., the peak around JD 2456250. 


B. Fe 11 LIGHT CURVES MEASURED BY 
INTEGRATION 

We present the Fe II light curves measured by the tradi¬ 
tional integration method for a comparison with that by the 
present fitting scheme in this Appendix. The flux of the Fe II 
bump on the red side of [O III] A5007 is calculated by inte¬ 
grating in 5115-5465 A the flux above the straight line set by 
two continuum windows 5085-5115 and 5465-5495 A (the 
three bands are in the rest frame). Figure[T0]shows these light 
curves. Comparing with those given by the fitting method 
(in Figures Q] and |5), all the light curves by integration show 
larger scatter. For four objects (Mrk 1044, Mrk 142, MCG 
+06-26-012, and IRAS F12397+3333), rough structures can 
be recognized in their light curves, and the derived Fe II time 
lags are consistent with those given by the fitting but with 
larger uncertainties. The remaining objects fail to have reli¬ 
able measurement of time lag, for their scattered light curves 
without clear structures. We also tried to measure the flux of 
the Fe II bump on the blue side of H/3 by choosing an inte¬ 
gration interval avoiding the broad He II line. Similarly, light 
curves with poor quality are obtained. 
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Figure 10. 
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Light curves of Fe 11 measured by simple integration (Fpe.intg)- All have larger scatter comparing with those obtained by fitting shown in Figures[4] 
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